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I.  IDENTIFICATION  AND  SIGNIFICANCE  OF  THE  PROBLEM  OR  OPPORTUNITY 

Speed  can  be  a  critical  advantage  lo  the  warfighter.  High-speed  propulsion  technologies  enable 
the  capability  to  defeat  a  time-critical  target,  to  escape  a  threat,  and  to  provide  rapid  reconnaissance  and 
responsive  access  to  space.  [1.2].  Rockets  provide  very  high  speeds,  but  with  limited  specific  impulse, 
which  is  the  duration  of  thrust  per  mass  of  propellant.  One  limiting  factor  in  rocket  performance  is  time- 
dependent  deformation  (creep)  in  the  throat  of  the  nozzle,  especially  in  solid  rockets  with  non-eroding 
throats.  On  the  other  hand,  hypersonic  aircraft  powered  by  scramjets  can  achieve  much  higher  specific 
impulse,  allowing  flight  durations  of  minutes  to  hours  at  speeds  in  excess  of  Mach  8.  However,  the 
performance  of  these  hypersonic  vehicles  is  limited  by  the  capability  of  the  vehicle  to  resist  the 
aerodynamic  heating,  particularly  on  the  leading  edges  of  the  wings  and  engine  inlet. 

Both  non-eroding  rocket  nozzles  and  sharp  leading  edges  for  hypersonic  vehicles  require 
materials  that  will  survive  service  conditions  over  2200*C  for  minutes  to  hours.  Meeting  these  needs  will 
require  significant  improvements  in  our  fundamental  understanding  of  the  response  of  materials  at  such 
high  temperatures,  especially  their  resistance  to  creep.  However,  conventional  methods  for  measuring 
creep  of  these  materials  have  severe  limitations  that  often  restrict  designers  to  extrapolating  properties 
from  lower-temperature  tests. 

Research  on  ultra-high  temperature  ceramic  (UHTC)  materials  has  been  performed  only 
sporadically  over  the  last  30  years,  resulting  in  “a  scarcity  of  fundamental  research  data  available  on  the 
mechanical  performance  and  manufacturability  of  UHTCs"  (3].  One  key  factor  limiting  the  maturation  of 
these  materials  is  our  inability  to  measure  their  performance  under  the  loads  and  at  the  extreme 
temperatures  to  which  they  are  exposed  in  service.  Our  team  has  demonstrated  a  new  method  for 
measuring  creep  at  up  to  2300*C  [4],  The  method  is  readily  extensible  to  over  3400*C  [5],  The  method 
has  been  demonstrated  at  stresses  exceeding  10  MPa,  with  higher  load  capability  under  separate 
development. 

Making  this  method  commercially  available  for  standardized  testing  of  UHTC’s  will  accelerate 
both  basic  research  and  the  development  of  new  materials  required  for  providing  greatly  enhanced 
capabilities  to  the  warfighter. 

2.  REVIEW  OF  PHASE  I  TECHNICAL  OBJECTIVES 

For  $100,000  and  9  months,  we  proposed  to  demonstrate  the  capability  of  our  method  to  measure 
creep  of  UHTC’s  at  temperatures  exceeding  2200®C.  using  small  inexpensive  samples  and  a  non-contact 
design.  Demonstration  tests  were  to  be  conducted  on  at  least  two  materials.  The  materials  selected  were 
pure  ZrB;  and  a  ZrBj  composite  with  25  volume  %  SiC,  The  samples  were  consolidated  from  10  jim 
ZrB:  powder  and  2  pm  SiC  powder  by  the  Naval  Surface  Warfare  Center.  Carderock  Division 
(NSWCCD),  as  described  by  Talmy,  er  a/.,  (6J, 

This  objective  was  divided  into  certain  research  questions;  Can  high-precision  UHTC  spheres  be 
processed  in  Electrostatic  Levitation  (ESL)  at  the  relevant  temperatures  and  for  the  required  duration? 
Can  ESL  provide  sufficient  stress  to  get  the  desired  strain  rates?  Can  the  demonstrated  methods  for  data 
collection,  reduction,  and  analysis  be  adapted  to  the  measurement  of  UHTC’s?  How  do  the  results 
compare  to  extrapolations  from  lower-temperature  conventional  tests? 

Proposed  tasks: 

1  Process  UHTC  materials  at  the  relevant  temperatures  in  Electrostatic  Levitation  for  extended 
periods. 

2,  Prepare  the  required  high-precision  spheres  from  UHTC  materials, 

3  Apply  sufficient  rotation  in  the  sample  to  cause  creep, 

4,  Measure  the  creep  deformation  to  sufficient  precision  and  accuracy  to  allow  analysis. 

5  Extract  the  constants  in  the  creep  constitutive  equation. 

6,  Compare  results  to  the  values  reported  in  the  literature  and  report. 
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3.  PHASE  I  WORK:  METHODS.  RESULTS,  AND  DISCUSSION 

3.1  Background  on  conventional  creep  methods 

Conventional  equipment  for  measuring  creep  uses  tensile  or/and  compressive  loads  applied  to 
samples  healed  by  an  electric  or  induction  furnace.  The  upper  end  of  a  specimen  is  gripped  by  a  fixture 
and  a  weight  is  attached  at  the  lower  end  to  apply  a  constant  load.  Fixtures  are  in  contact  with  the 
specimen,  therefore,  they  must  be  made  of  the  materials  which  have  equal  or  better  chemical  and 
mechanical  stability  than  the  specimen  at  high  temperatures.  At  high  temperatures,  the  specimen  becomes 
more  reactive,  such  that  it  may  be  incompatible  with  the  testing  apparatus.  Most  of  the  conventional 
methods  for  the  measurements  of  creep  are  limited  to  below  -1. 700  *C.  For  testing  in  the  range  of 
1400*C  to  1700*C.  a  thermal  gradient  design  is  required  to  cool  the  grips.  This  arrangement  requires 
extremely  long  specimens  which  are  expensive  and  time  consuming  to  fabricate  and  finish.  Notable 
exceptions  are  measurements  of  tungsten  and  its  alloys  at  temperatures  up  to  2800  "C,  using  electric 
resistance  healing  in  the  specimen  and  large  thermal  gradients.  However,  these  methods  are  not  applicable 
to  nonconductors  or  to  materials  which  cannot  be  formed  into  a  wire. 

Creep  is  the  time-dependent  deformation  of  a  material  under  constant  stress,  usually  at  high 
temperature.  The  strain  rate  s  has  been  observed  to  follow  the  relation  e  =  Ca'‘ f{r^  where  a  is  the 
stress,  n  is  the  stress  exponent.  f{T)  is  a  function  of  temperature,  and  C  is  a  constant.  The  stress 
exponent  n  is  determined  from  multiple  measurements  at  different  stresses  at  the  same  temperature,  with 
a  typical  accuracy  of  many  tens  of  percent  in  n  (Fig,  1.) 

Because  of  the  limitations  of  conventional  testing  methods,  high-temperature  creep  behavior  is 
often  extrapolated  from  low-temperature  measurements.  This  extrapolation  is  often  inaccurate,  especially 
for  multi-phase  materials,  A  new  non-contact  method  for  measuring  creep  of  ultra-high-temperature 
(>2,000  *C)  materials  is  strongly  demanded.  This  method  must  be  to  both  metallic  and  nonmetallic 
materials,  for  the  future  technologies  of  high  performance  rockets,  hypersonic  aircraft,  and  high  efficiency 
energy  conversion  systems, 
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Figure  I  Strain  rate  vs.  stress  for  conventional  creep  measurements, 
Nbat  1,985  "C  under  1.0, 1.5,  and  2.0  MPa. 
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3.2  Non-contact  creep  method 

Our  team  has  developed  a  novel  approach  to  measuring  creep  at  extremely  high  temperatures 
using  electrostatic  levitation.  This  method  has  been  demonstrated  to  measure  creep  of  niobium  up  to 
2300*C  [4],  while  ESL  has  melted  tungsten  (3400*C)  [5],  The  small  samples  used  in  this  method.  2-3 
mm  diameter  spheres,  allow  a  significant  advantage  in  material  utilization  and  cost  of  fabricating  samples. 

The  basis  of  this  method  is  that  a  rotating  sphere  experiences  a  shear  stress  r  that  is  determined 
by  its  angular  velocity  o,  radius  R ,  and  density  p .  The  maximum  shear  stress  is  =  O.lWorR^ p  . 
The  sample  then  creeps  under  the  applied  load.  ForZrBj  rotating  at  4000  Hz  (240.000  RPM).  the  shear 
stress  developed  at  the  center  is  approximately  1  MPa, 

Achieving  such  high  rotation  rates  is  challenging,  but  not  impossible.  Since  the  ESL  samples  are 
heated  by  a  laser,  we  have  used  photon  pressure  to  accelerate  the  samples.  Each  absorbed  photon 
transfers  its  momentum  to  the  sample.  The  heating  laser  is  aligned  off  the  axis  of  the  sample,  so  that  this 
momentum  transfer  results  in  a  net  torque.  For  2,2  mm  diameter  ZrBi  spheres  at  2200®C,  using  the 
properties  in  references  [7,8]  gives  an  estimate  of  about  45  minutes  to  reach  the  4000  Hz  required  for  1 
MPa  stress,  10  MPa  is  predicted  to  require  12.650  Hz  and  just  over  2  hours  of  acceleration.  For 
comparison,  our  experiments  achieving  just  less  than  1  MPa  in  niobium  required  about  3000  Hz  rotation. 

Since  summer  2009,  most  experiments  have  been  conducted  using  an  induction  motor  to  cause 
the  rotation  of  the  samples  instead  of  photon  pressure.  The  induction  motor  decouples  load  from 
temperature,  allowing  higher  stresses  at  lower  temperatures,  but  is  limited  to  conducting  samples. 

As  the  sample  deforms,  its  surface  shape  is  captured  by  high-speed  video  (Figs,  2-3),  The  shape 
at  each  lime  is  measured  by  a  high-precision  machine  vision  algorithm  developed  at  the  University  of 
Massachusetts  [9].  This  method  can  determine  the  surface  shape  to  approximately  170  parts  per  million 
(ppm),  or  approximately  190  nm  [10].  The  defonnations  are  compared  to  finite  element  models  to 
determine  the  constitutive  constants  in  the  creep  relation.  Furthermore,  the  shape  of  the  surface  is 
uniquely  determined  by  the  stress  exponent.  With  our  method,  the  stress  exponent  can  be  determined  to 
an  accuracy  of  a  few  percent  and  with  a  single  test.  This  accuracy  compares  favorably  with  the  many  lens 
of  percent  typical  of  conventional  methods,  which  also  require  many  tests  at  different  stress  levels  to 
determine  the  stress  exponent. 
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Figure  2  Captured  images  of  deformed  sample.  Equatorial  radius  increases  and  polar  radius 

decreases  as  the  deformatiou  proceeds. 
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Figure  3  Plot  of  tbe  change  both  polar  and  equatorial  radii  of  tbc  sample.  In  tbc  beginning,  the 
equatorial  radius  deercases  due  lo  the  mass  evaporation. 

Our  non-contoct  method  was  validated  in  collaboration  with  the  University  of  Tennessee  for 
niobium  at  1985  ‘’C.  Through  the  heroic  efforts  of  our  contractor.  Pittsburgh  Materials  Testing.  Inc,, 
conventional  tests  were  completed  at  this  temperature.  Conventional  tests  at  2000®C  were  unsuccessful 
due  to  a  reaction  between  the  sample  and  the  grips.  Comparison  with  conventional  experiments  showed 
agreement  within  the  uncertainty  of  the  conventional  measurements. 


Tabic  1  Stress  exponent  of  Nb  at  1985‘'C 


ESL.  UMass/NASA,  2006 

2.52  ±0.02 

Conventional  Test,  UT,  2006 

2.4  ±1 

Kiessig  &  Essniann,  (11]  1985 

2.476  ±? 

Frost  &  Ashby.  |12|  1982 

4.4  ±? 

33  Organization  and  Project  Management 

The  project  is  a  collaboration  among  Dr,  Robert  W,  Hyers  (RHA,  PI),  Dr.  Jan  R,  Rogers  {NASA,  co-1), 
and  Prof  Mario  Rotea  (UMass,  co-1).  RHA  managed  the  overall  project  and  is  responsible  for  delivery  of 
the  deliverables  proposed,  NASA  provided  ESL  processing  and  development  of  processing  techniques  in 
accordance  with  the  subaward.  UMass  provided  data  reduction  and  analysis  and  finite  element  modeling 
in  accordance  with  the  sub  award.  Other  personnel  involved  in  the  project  were  project  manager  Liz 
Hyers.  RHA;  graduate  students  Stacy  Cane  pari  and  Xiao  Ye.  UMass;  and  ESL  staff  Trudy  Allen  and 
Glenn  Fountain  at  NASA 
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3.4  Task  1:  Process  UHTC  materials  at  the  relevant  temperatures  in  Electrostatic  Levitation  for 
extended  periods. 

Preliminary  work  on  this  task  had  been  performed  in  September  2tX)6  at  the  NASA  MSFC  ESL 
by  Dr,  Jan  Rogers.  Tom  Rathz,  and  Trudy  Allen,  Pure  ZrB;  material  was  donated  by  Dr,  Eric  Wuchina  of 
the  NSWCCD.  The  NSWCCD  provided  more  ZrB-  and  ZrBj-SiC  composite  material  for  the  current 
effort. 

In  2006,  a  hand-filed  spheroid  of  this  material  was  processed  in  ESL  at  up  to  2500®C,  healed  by  a 
carbon  dioxide  laser  (10.6  micron  wavelength).  The  temperature  was  measured  by  pyrometry,  and  the 
relevant  spectral  emissivities  were  recorded  by  the  ESL  team.  The  sample  showed  an  insignificant  mass 
loss  of  approximately  0.5%.  Residual  gas  analysis  detected  strong  peaks  of  Helium,  consistent  with  the 
processing  of  the  sample  material.  No  special  handling  was  required.  The  samples  were  held  over 
I950*C  fora  total  of  approximately  one  hour.  Additional  processing  in  the  spring  of  2009  confirmed  the 
processability  of  low-precision  samples  of  this  material. 

Task  1  was  completed  successftilly;  processing  of  ZrBj  in  ESL  at  temperatures  over2200*C  has 
been  demonstrated. 

3.5  Task  2:  Prepare  the  required  bigh-prccision  spheres  from  UHTC  materials. 

RHA  suhconli'acled  the  task  of  grinding  the  spherical  samples  to  of  Industrial  Tectonics,  Inc, 

(ITl)  of  Dexter.  MI  (www.itibalLcom),  The  NASA  ESL  and  UMass  groups  have  employed  Ill’s  services 
for  fabrication  of  high-precision  spheres  for  many  years.  The  4  weeks  estimated  in  the  proposal  for 
sample  fabrication  in  the  proposal  proved  too  optimistic;  in  contrast  to  our  prior  experience,  negotiation 
of  the  subcontract  alone  took  that  long.  ITI  experienced  several  issues  with  sample  fobrication,  related  to 
roundness  and  surface  finish  of  the  spheres.  In  the  end,  ITl  was  able  to  manufacture  ball  grade  25  spheres, 
with  a  sphericity  of  25  microinches  (-635  nm),  but  not  until  after  the  period  of  performance  of  this  Phase 
I  project  had  expired.  Experience  now  indicates  that  six  months  should  be  allowed  for  ^bricalion  of 
samples  from  materials  like  ZrB;,  RHA  has  also  identified  a  potential  second  source,  PremaTech 
Advanced  Ceramics  in  Worcester.  MA,  who  may  also  be  able  to  manufacture  high-precision  spheres  for 
future  experiments. 

Task  2  was  completed  successfully:  high-precision  samples  were  successfully  fabricated. 

3.6  Task  3:  Apply  sutllcicnt  rotation  in  the  sample  to  cause  creep. 

The  preliminary  estimates  given  in  the  proposal  in  February  2008  were  based  on  a  stress  of  1 
MPa.  corresponding  to  a  rotation  rate  just  under  4000  Hz,  The  flexural  creep  measurements  of  Talmy.  et 
al.,  [6],  published  in  May  2008,  show  that  much  higher  stresses  are  required  to  get  the  required 
deformation  within  a  reasonable  experiment  duration.  Finite  element  calculations  by  our  research 
institution.  UMass.  confirm  that  a  rotation  rate  closer  to  32,000  Hz  is  required  to  provide  the  100  MPa  of 
stress  required  to  get  strain  rates  of  the  order  of  I  O'Vsec  to  allow  the  sample  to  deform  -9%  within  a  l- 
day  test  (13]  (attached  as  Appendix).  At  this  rotation  rate,  sufficient  deformation  is  predicted  for 
temperatures  exceeding  1800"C  for  pure  ZrB;,  and  15()0®C  for  ZrB;  +  25  vol%  SiC,  The  lower  density  of 
the  composite  does  cause  a  lower  stress  at  the  same  rotation  rate,  but  this  lower  stress  is  more  than 
compensated  for  by  the  faster  creep  of  the  composite. 

Meanwhile,  a  magnetic  device  for  inducing  sample  rotation  was  developed  under  separate 
support  through  NASA’s  Innovative  Partnerships  Program,  This  new  method  uses  a  rotating  magnetic 
field  to  induce  a  torque  in  electrically  conducting  samples.  It  was  demonstrated  successfully  on  low- 
precision  ZrB;  spheres  up  to  4000  Hz  in  July  2009,  limited  by  problems  related  to  the  shape  of  these  low- 
precision  samples.  The  magnetic  rotation  device  has  since  been  demonstrated  to  cause  rotation  rates  over 
13,000  Hz  on  high-precision  spheres  of  other  materials.  The  hardware  was  designed  to  induce  rotation  at 
rates  up  to  and  exceeding  the  32,000  Hz  indicated  by  the  finite  element  simulations. 

Coupled  with  the  latest  measurements,  simulations  show  that  rotation  rates  near  32.000  Hz  are 
required  to  a  cause  sufficient  creep  rate  to  complete  the  experiment  in  one  day.  The  electromagnetic 
induction  motor  was  designed  to  provide  rotation  rates  this  high.  Experiments  show  that  the 
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electromagnetic  induction  motor  couples  well  to  the  ZrBi  samples.  These  factors  indicate  that  creep 
measurements  using  this  method  are  feasible  in  ESL,  even  though  the  integrated  demonstration  of  creep 
originally  planned  was  delayed  past  the  end  of  the  period  of  performance  by  the  problems  with  sample 
manufocture, 

3.7  Task  4:  Measure  the  ereep  deformation  with  sufficient  precision  and  aecuracy. 

In  measuring  the  creep  deformation,  a  few  requirements  must  be  met.  First,  the  sample  must 
deform  axisymmetrically.  Accordingly,  the  samples  must  start  axisymmetric,  so  high-precision  spheres 
are  used  for  these  measurements.  Achieving  symmetrical  deformation  also  requires  that  the  samples  be 
either  single  crystals  or  fine  polycrystals.  The  materials  received  from  the  NSWCCD  are  of  the  latter 
category;  with  approximately  200  grains  spanning  the  sample  diameter,  they  are  expected  to  deform 
isotropically. 

Second,  the  video  technique  must  be  able  to  determine  the  surhice  shape  very  precisely.  We 
estimate  the  surface  was  determined  to  approximately  170  parts  per  million  (ppm),  or  approximately  190 
nm.  Such  accuracy  is  necessary  to  allow  determination  of  the  stress  exponent  to  an  accuracy  of  a  few 
percent  from  a  single  test.  Besides  axisymmetry,  the  sample  must  also  have  a  smooth  surface  and  provide 
sufficient  contrast  for  imaging  the  surface  with  the  high-speed  camera.  The  smooth  surface  is  provided 
by  the  grinding  operations  by  ITl,  At  the  target  measurement  temperatures,  the  samples  are  sufficiently 
self-illuminated  to  provide  good  contrast  as  well.  The  video  imagery  will  be  acquired  with  NASA  ESL’s 
Phantom  V7  high-speed  digital  camera.  The  use  of  this  camera  is  included  with  the  use  of  their  facility. 

Successfiil  completion  of  this  task  required  high-precision  spheres  that  were  not  available  within 
the  period  of  performance  of  this  project. 

3.8  Task  S:  Extract  the  constants  in  the  creep  constitutive  equation. 

This  task  involves  running  finite  element  simulations  of  the  deformation  of  hypothetical  samples 
under  the  load  history  experienced  by  the  actual  sample,  and  using  those  results  to  determine  the 
parameters  in  the  constitutive  equation  for  creep  that  best  match  the  response  of  the  actual  sample.  This 
task  was  intended  to  be  applied  to  the  results  of  Task  4.  Lacking  those  results,  the  simulations  were 
instead  performed  using  the  same  finite  element  models  to  predict  the  performance  under  different  ideal 
experimental  conditions.  These  results  are  summarized  in  section  3.6,  above,  and  reported  in  detail  in 

[13]. 

The  feasibility  of  this  task  had  been  established  on  Nb  samples  prior  to  the  Phase  1  work,  and  was 
demonstrated  in  simulation  during  the  Phase  I  project. 

3.9  Task  6:  Compare  results  to  the  values  reported  in  the  literature  and  report. 

The  results  of  the  simulations  were  presented  at  the  AFOSR  Workshop  on  Materials  for  Extreme 
Environments.  August  3-5,  2009,  and  submitted  to  the  Journal  of  the  European  Ceramic  Society  [13], 
These  simulations  show  that  we  can  expect  to  compare  directly  to  the  highest  temperature  measured  by 
Talmy,  et  al ,  [6]  for  the  composite,  but  that  we  can  only  compare  to  extrapolations  of  their  results  for  the 
pure  ZrB:,  The  preliminary  experiments  with  low-precision  spheres  were  also  presented  at  the  workshop, 

4  SUMMARY  AND  CONCLUSION 

A  9-month  study  of  the  feasibility  of  applying  the  non-contact  technique  for  creep  measurement 
to  ZrB;  and  ZrB;  +  SiC  composites  was  carried  out.  Processing  of  low-precision  samples  shows  that 
these  materials  may  be  processed  in  ESL  at  temperatures  over  2000"C,  and  up  to  2500®C  for  the  pure 
material.  The  composite  shows  a  eutectic  at  2270"C,  and  so  should  not  be  processed  too  close  to  this 
temperature.  The  materials  couple  well  with  the  new  electromagnetic  rotator,  allowing  much  higher 
stresses  as  well  as  expanding  the  time  under  load  that  may  be  accomplished  in  the  experiments.  Finite 
element  simulations  were  performed  with  parameters  extrapolated  from  reported  measurements  of  the 
creep  of  ZrB^  and  ZrBi  +  SiC  composites.  These  simulations  indicate  that  the  experiments  are  feasible 
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within  the  designed  performance  of  the  electromagnetic  rotator  within  the  allotted  experiment  duration. 
High-precision  samples  of  ball  grade  25  were  successfully  fabricated  of  both  materials,  ZrB^  and  ZrB;  + 
25  vol%  SiC.  These  factors  combine  to  provide  a  strongly  positive  indication  of  the  feasibility  of 
applying  this  method  to  these  ultra-high-temperature  materials. 

A  series  of  delays,  including  a  time  to  manufacture  the  samples  that  was  5  months  longer  than 
projected,  meant  that  the  final  integrated  test  with  high-precision  spheres  could  not  be  completed  within 
the  period  of  performance  of  the  contract. 
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APPENDIX:  MANUSCRIPT  DESCRIBING  SIMULATION  METHODOLOGY  AND  RESULTS 

The  following  manuscript  is  reference  [13],  which  was  submitted  in  September  2009  for  review 
by  the  Journal  of  the  European  Ceramic  Society, 
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Computational  Methods  for  the  Analysis  of 
Non-contact  Creep  Deformation 
in  ZrBi-SiC  Composites 

Xiao  Ye  and  Robert  W.  Hye^s^ 

Deparlmeni  of  Mechanical  and  Industrial  Engineering,  University  of  Massachuseils, 
160  Governors  Dr.  Amherst,  Massachusetts  OHiOS.  USA 

Abstract: 

A  need  for  higher  service  temperatures  is  driving  development  of  new,  higher- 
temperature  materials  that  are  resistant  to  creep.  However,  conventional  methods  of 
measuring  creep  become  increasingly  difficult  over  about  1700®C.  A  non-contact 
method  with  the  capability  of  measurements  at  much  higher  temperature  has  been 
demonstrated  on  niobium  using  centrifugal  loading  of  a  spherical  sample.  Recent  efforts 
have  been  made  to  extend  this  method  to  lower  temperatures  and  higher  stresses.  Using 
material  properties  from  the  literature,  we  performed  finite  element  simulations  to 
determine  the  range  of  experimental  parameters  over  which  non-contact  measurements  of 
creep  can  be  readily  finished  within  a  reasonable  time  duration  for  ZrB^  and  ZrB^  +  25 
vol%  SiC.  Results  from  FEA  model  shows  that  the  experiments  are  feasible  at  an  angular 
velocity  of  32,000  rps  and  a  temperature  as  low  as  1 800  for  ZrB2  and  1500  for 

ZrBs  +25  vol%SiC, 


'  CorrcsiHinding  author.  Tel  +1-413-545-2253;  fax +1  -4 13-545- 1027.  E-mail  address. 
hycrs^ecs.umass.cdu. 
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1.  Introduction 

Modem  applications  such  as  rocket  nozzles,  leading  edges  of  hypersonic  airplanes, 
and  next  generation  turbine  blades  place  demanding  requirements  on  ultra-high 
temperature  ceramics.  Besides  requirements  for  general  mechanical  strength  and 
oxidation  resistance,  at  high  temperature,  creep  resistance  is  one  of  the  most  crucial 
criteria.  UHTCs,  including  transition  metal  carbides,  oxides  and  borides  are  of  current 
interest  for  these  applications.  Zirconium  diboride  (ZrBj)  is  characterized  by  a  high 
melting  point  (3040®C),  high  hardness  {22  GPa)  and  good  corrosion  resistance  due  to  its 
high  percentage  of  covalent  bonding,'  However,  at  temperature  greater  than  650  ®C,  its 
resistance  to  chemical  attack  by  oxidation  is  not  satisfactory  to  accommodate 
atmospheric  applications.  In  order  to  improve  its  oxidation  resistance  at  elevated 
temperature,  different  additions  have  been  introduced.  Among  them,  25  vol%  SiC  is 
currently  the  most  promising  for  applications  in  oxidizing  environments.^ 

/  I  Fundamental  creep  behavior 

Creep  is  defined  as  a  time-independent  deformation  that  usually  should  be  included 
into  consideration  when  using  a  material  in  conditions  like  this:  temperature  higher  than 
half  of  its  absolute  melting  point  and  subjecting  to  a  stress  load  lower  than  its  yielding 
point.  The  rate  of  creep  is  affected  by  various  factors,  such  as  material  properties, 
exposure  temperature,  exposure  time  and  applied  stress.  Despite  the  differences  between 
various  creep  mechanisms,  a  complete  creep  behavior  generally  can  be  divided  into  three 
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stages  based  on  their  different  creep  rate  characteristics.  In  stage  1 ,  (transient  creep),  the 
creep  rate  decreases  until  a  minimum  steady-state  is  reached,  which  is  stage  2. 
Transmission-microscopic  examination  shows  that  at  the  beginning  stage  material 
structure  evolves  with  a  significant  strain  accumulation,  which  corresponds  to  a 
dramatically  dislocation  density  rise.  And  then  during  the  following  stage,  speed  of  new 
dislocation  generation  and  annihilation  among  dislocation  balance  to  each  other  to  reach 
a  steady  state  macroscop ically.  In  stage  3,  creep  rate  rise  again  until  fracture  happens. 

1.2  Non-contaci  measurement  of  creep 

Since  creep  resistance  is  an  important  criterion  in  demanding  applications,  accurate 
measurement  of  creep  behavior  of  candidate  materials  will  be  required.  At  high 
temperature,  materials  are  highly  reactive.  Reaction  with  the  test  fixtures  generally 
excludes  conventional  method  for  measuring  creep  above  about  1700  ®C.  A  group  in 
UMASS  has  been  working  on  non-contact  measurement  of  creep  behavior  since  2004, 
Lee^^ii/.^  used  Electrostatic  Levitation  (ESL)  to  levitate  spherical  Nb  specimen  under 
high  vacuum  ( ^  10’’  torr)  at  NASA  Marshall  Space  Flight  Center  (MSFC)  in  Huntsville, 
AL.  This  apparatus  has  a  demonstrated  capability  to  process  samples  at  extremely  high 
temperatures,  having  been  used  to  melt  tungsten  (3410 

Specimens  were  stabilized  via  three  pairs  of  orthogonal  electrodes  controlled  by  PID 
control-loop,  A  200  W  YAG  laser  beam  was  cast  on  the  sphere  to  heat  and  rotate  the 
sample.  The  shape  of  the  spherical  sample  was  captured  with  a  digital  camera  and  the 
captured  images  were  analyzed  to  determine  the  amount  of  deformation  of  the  sample. 
Meanwhile  an  ANSYS  model  was  established  to  simulate  creep  behavior  of  the  sphere 
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with  different  stress  exponent  at  the  same  conditions  as  real  experiment.  By  comparing 
the  deformed  shape  of  the  sample  to  the  AN  SYS  model,  the  stress  exponent  of  the 
specimen  can  be  determined,  which  further  helps  in  determine  the  creep  mechanism.  The 
stress  exponent  measured  in  this  way  agrees  with  conventional  tensile  measurement. 

In  this  paper,  based  on  the  same  idea,  creep  data  for  ZrB^  and  ZrBj  +  25  vol%  SiC 
was  collected  by  screening  available  literature  and  input  into  the  model.  Ideal  simulation 
results  will  help  evaluate  the  feasibility  of  using  this  technique  for  UHTC’s  and  to  plan 
the  experimental  matrix, 

/.i  Experimentally  determined  creep  behavior  in  UHTCs 

I,  G.  Talmy  et  al?  characterized  flexural  creep  of  ZrB^  +  ()“50  vol%  SiC  as  a  function 
of  temperature  (1 200-1 500  ®C},  stress  (30-180  MPa)  and  SiC  particle  size  (2  and  \Q  pm) 
in  ambient  air  via  four-point  flexural  creep  test.  Specimens  were  synthesized  by  hot 
pressing  commercially  available  ZrB^  and  SiC  power  at  2100  ®C  for  1  h  with  20  MPa 
pressure  in  He  atmosphere.  The  creep  rate  increased  with  increasing  SiC  content, 
temperature,  and  stress  and  with  decreasing  SiC  particle  size,  especially,  at  temperatures 
above  1300  ®C.  The  activation  energy  of  creep  turned  out  to  be  a  function  of  SiC  content 
and  increased  linearly  with  it  for  ceramics  containing  0-50  voI%  lpmS\C.  At  1400  ®C, 
the  reported  stress  exponent  was  1.1  for  Pure  ZrBj  and  0,7  for  ZrB2/25  vol%  2  ^mSiC 
respectively,  which  indicated  for  both  compositions  diffusional  creep  had  a  great 
contribution  to  total  creep  deformation  at  these  temperatures.  Cracking  and  grain  shifting 
were  observed  on  the  tensile  side  of  25  vol%  SiC  specimen.  Also,  the  presence  of  stress, 
both  compressive  and  tensile,  promoted  oxidation. 
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Melendez-Martinez  et  a  A*  evaluated  creep  behavior  of  pure  ZrBz  and  ZrBj  +  4  wt% 
Ni  at  lemperarures  between  1400  and  1600  ®C  and  at  stresses  ranging  from  47.0  to  472.3 
MPa  for  pure  ZrB^  and  10  to  63,5  MPa  for  Ni-doped  ZrBj,  The  compressive  creep  test 
was  performed  in  controlled  argon  atmosphere  via  a  prototype  dead- weight  load  device 
on  specimens  which  were  hot  pressed  at  1900  ®C  (pure  ZrB^)  and  1850  ®C  (Ni-doped 
ZrBj)  for  30  min  with  applied  pressure  of  30 MPa,  The  following  mechanical  properties 
were  measured:  elastic  modulus  (E)  by  the  resonance  frequency  method,  flexural  strength 
(S)  by  four-point  flexure  test,  microhardness  (Hvl  .0)  by  a  Vickers  indenter.  Pure  ZrBj 
showed  only  little  creep  deformation  at  temperature  below  1400  and  stress  below 
298 MPa.  The  reported  stress  exponent  was  1,7  at  1500^^0  and  0,6  at  1600®C  below 
220 MPa.  The  Ni-doped  ZrBj  failed  catastrophically  for  stress  high  than  2 5 MPa  and 
shows  a  ductile  behavior  only  at  lower  stress.  The  stress  exponent  was  1 .5  at  1500  ®C  and 
stress  between  10  and  20  MPa,  This  behavior  may  be  attributed  to  the  presence  of  Ni-rich 
grain  boundary  phases  at  triple  points  of  the  grain  structure. 

2.  Procedure 

A  sphere  model  with  radius  1, 1mm  was  established  in  ANSYS.  Creep  behavior  of 
two  different  materials  was  simulated  with  ANSYS:  Pure  ZrBj  and  ZrB^  +  25  vol%  SiC. 

An  implicit  creep  equation  was  chose  to  model  the  creep  behavior  of  the  sphere. 
Values  for  parameters  used  in  this  equation  were  collected  by  surveying  the  literature. 
Some  of  them  were  well  known  or  could  be  easily  extrapolated  fi'om  available  data,  such 
as  density  ( p\  Young’s  modulus  (E),  Poisson’s  ratio  (t; ),  others  were  extrapolated  from 
existing  data.  Since  the  purpose  of  this  model  is  to  evaluate  feasible  range  for 
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experiments,  and  the  uncertainty  involved  in  these  extrapolations  is  acceptable  in  the 
absence  of  previous  measurements.  These  material  properties  were  input  into  the  finite 
element  model  of  the  rotating  sample.  The  model  was  run  for  various  conditions  to 
determine  the  range  of  experimental  parameters  that  is  expected  to  provide  sufficient 
creep  deformation  within  the  duration  of  experiment. 

2J  Mec  hanic  al  properties 

Recently,  ZrB^  and  ZrB2-based  composites  have  been  densifled  by  various  methods, 
including  hot  pressing,  reactive  hot  pressing,  plasma  sintering  and  pressureless  sintering. 
Slight  changes  in  processing  conditions,  such  as  temperature,  thermal  history,  hold  time, 
pressure,  or  species  and  quantities  of  additives  will  all  have  reflections  in  product's 
mechanical  properties. 

Chamberlain  ei  al^  densified  ZrB^  and  ZrBj  ceramics  containing  10, 20  and  30  vol% 
Sic  particulates  from  commercially  available  powders  by  hot  pressing  and  measured 
density,  four-point  bend  strength,  fracture  toughness,  elastic  modulus  and  hardness  at 
room  temperature.  The  data  for  ZrBj  +  25  vol%  SiC,  can  be  interpolated  by  taking 
arithmetic  mean  of  the  same  properties  of  materials  with  close  composition.  Table  1 
shows  the  data  for  ZrBj  and  ZrB^  +  25  vol%  SiC  at  room  temperature. 

At  high  temperature,  it  is  hard  to  measure  properties  like  Young’s  modulus.  There  are 
only  limited  data  of  mechanical  properties  concerning  transition  metal  borides.  Opeka  et 
at}  measured  mechanical  properties  for  HfB2,  HfCo^s,  HfCos?,  and  HfNo^i  ceramics. 
They  reported  that  all  these  ceramics  will  experience  a  dramatic  modulus  loss  beginning 
as  low  as  1090  HfBj  lost  about  75%  of  its  room  temperature  modulus  at  1640 
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Accordingly,  in  order  to  account  for  this  modulus  loss,  a  Young's  modulus  of  120  GPa 
was  used  for  pure  ZrBj  and  130  GPa  for  ZrB^  +  25  voI%  SiC  for  all  the  temperatures 
considered  in  our  simulation.  Since  Young's  modulus  would  have  a  considerable  impact 
on  only  the  elastic  deformation  of  specimen  sphere,  which  is  much  smaller  than  the  creep 
deformation  (plastic  deformation),  the  accuracy  of  this  estimation  should  be  sufTicieot, 
Another  mechanical  property  required  for  ANSYS  model  is  Poisson's  ratio.  In 
contrast  to  other  properties,  it  is  insensitive  to  additives  as  well  as  to  porosity  caused  by 
different  processing  conditions.  It  remains  almost  constant  for  various  ZrBj-based 
composites.®*'®  Guo  et  at.  measured  mechanical  properties  for  ten  different  series  of  hot 
pressed  ZrB2-MoSi2“SiC,  ranging  from  10  to  40  vol%  for  MoSij,  and  5  to  20  voI%  for 
Sic,  at  1 800  ®C  for  30  min  under  a  pressure  of  20MPa  in  vacuum.  Poisson’s  ratio  barely 
changes  for  these  different  materials.  Table  2  shows  Poisson’s  ratio  for  pure  ZrBj  and 
ZrBs  +  25  vol%  SiC  we  used  in  all  simulations.  Small  changes  due  to  differences  from 
the  reference  temperature  will  be  neglected. 

2.2  Creep  model  and  determination  of  its  parameters 

As  a  powerful  and  widely  used  Finite  Element  Analysis  (FEA)  software  package, 
ANSYS  was  used  to  simulate  creep  behavior  of  monolithic  ZrB2  and  ZrB;  +  25  voI%  SiC. 
ANSYS  provides  as  much  as  13  featured  implicit  creep  models  to  simulate  primary  and 
secondary  creep  behavior.  The  Norton  model  (eq.  1}  was  chosen  for  these  simulations  due 
to  its  programming-friendly  concise  expression: 

^  =  Ccr"i‘^)  (1) 
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where  ^  is  the  strain  rate,  coefficient  C  is  a  material  constant,  n  is  the  stress  exponent,  Q 
is  the  activation  energy,  and  R  is  the  universal  gas  constant,  which  equals 

f  mol ■  K .  The  stress  exponent  is  determined  by  the  creep  mechanism,  which  is 
itself  dependent  on  temperature.  For  diffusional  creep  usually  dominant  at  low  stress, 
such  as  Coble  creep  and  Nabaito-H erring  creep,  n  is  close  to  1,  while  for  power  law 
creep,  which  usually  happens  at  high  stress,  n  is  greater  than  1.  The  thermal  activation 
energy,  in  contrast,  is  insensitive  to  creep  mechanism  but  depends  mainly  on  material’s 
composition.  The  thermally  activated  nature  of  creep  implies  a  strong  dependence  on 
temperature  of  the  strain  rate, 

1,  G.  Talmy  et  al?  reported  that  stress  exponent  was  about  1  for  ZrB^  containing  0-25 
vol%  SiC  around  1400  ®C,  which  indicated  that  diffusional  creep  was  the  dominant  creep 
mechanism.  The  activation  energy  didn't  change  within  test  temperature  ranging  from 
1200  to  1500  ®C  for  ZrBj -based  composites  of  varying  SiC  content:  1 30/t/  /  for 
pure  ZrBi  and  21(ikJ  I  mole  tax  ZrB2.+  25  vol%  SiC.  The  leading  coefficient  C  can  be 
calculated  using  equation  (1),  which  yields  3. 965 x  10*  /  sec  for  pure  ZrB^  and 
2.084X  10  ^  /  sec  for  ZrB2.+  25  vol%  SiC. 

Lee  et  at}  used  the  same  creep  model  for  simulating  niobium  (Nb)  sphere's  creep 
behavior.  Angular  velocity  was  set  to  be  20106  rad/s  (3200  rps)  to  apply  a  shear  stress 
load  about  I  MPa  at  sphere  center.  They  reported  a  3300s  predicted  duration  of 
experiment  to  reach  the  stop  criterion  of  9%  equatorial  creep  deformation,  which  was  the 
same  as  we  used  here.  The  creep  rate  was  at  the  order  of  10*^5  ’ . 

The  maximum  shear  stress  in  the  sample  varies  with  sample  size  and  rotation  rate 
according. 
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r  =  0.2nft.VV  (2) 

where  <y  is  angular  velociry  (rad/s},  r  is  radius  of  the  sphere,  p  is  material  density. 

Achieving  stresses  similar  to  those  reported  by  Talmy  et  al.^  requires  an  angular 
velocity  of  201060  rad/s  (32,000rps). 

However,  the  creep  rate  reported  by  Talmy  ei  o/.  at  1400  and  lOOMPa  pressure  is 
only  10  “.v"'  for  pure  ZrBj  and  10  ’.v  '  for  ZrB2  +  25  vol%  SiC,  respectively,  both  of 
which  are  too  small  to  reach  a  satisfactory  strain  of  9%  equatorial  deformation  within  a 
few  hours.  Increasing  temperature  is  necessary  to  have  a  creep  rate  comparable  to  Lee's. 
Extrapolation  of  Talmy 's  data  gives  1800  ®C  for  pure  ZrBj  and  1 500  ®C  for  ZrBj  +  25 
vol%  Sic.  The  corresponding  stress  exponent  and  resulting  coefficient  C  from  equation 
i  1 )  are  given  in  table  3. 

3.  Result  and  discussion 

For  every  1000  simulation  seconds  before  94%  of  desired  strain  and  later  on  10 
simulation  seconds,  angular  velocity  was  updated  to  account  for  changes  in  polar  moment 
of  inertia  which  resulted  from  shape  change.  The  simulation  runs  about  3  hours  for  pure 
ZrB;  and  5  hours  for  ZrBj  +  25  vol%  SiC  in  an  Intel  core  2  computer. 

For  pure  ZrBj,  the  model  shows  that  the  equatorial  and  polar  radii  are  L199x  10 'm 
and  9,287  X 1  O'*  m,  respectively,  which  corresponds  to  9%  strain  at  the  equator  after  an 
elapsed  time  of  5102  seconds.  The  deformation  in  thex  direction  of  the  XZ  plane  and 
von  Mises  stress  contour  are  plotted  in  Figures  1  and  2,  respectively. 
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For  ZtBz  +  25  vol%  SiC,  model  showed  lhal  the  equatorial  and  polar  radii  are 
1.199x10  ‘m  and9,164xl0‘^m,  respectively,  at  run  time  12272  sec.  X  direction 
deformation  &  von  Mises  stress  contour  are  showed  in  Figures  3  and  4. 


The  simulations  show  that  completing  experiments  with  an  8-hour  shift  is  feasible,  if 
only  at  very  high  angular  velocity  around  32,000  rps.  Although  the  method  has  been 
demonstrated  only  at  speeds  up  to  1 1,000  rps,  work  is  underway  at  UMASS  and  NASA 
to  extend  the  experimental  range  to  speeds  over  30,000  rps.  The  results  of  the  simulations 
here  show  the  importance  of  this  capability  for  higher  speeds  to  enable  the  high  loads 
needed  for  testing  creep  in  UHTC’s  in  a  short-duration  experiment.  With  the  extended 
range  of  rotational  speed  and  stress,  the  range  of  applicability  of  the  new  method  overlaps 
with  conventional  methods,  while  still  providing  the  capability  of  measurements  at  much 
higher  temperatures  than  conventional  methods. 


4.  Summary 

A  sphere  model  was  established  in  ANSYS.  Rotating  at  an  angular  velocity  of 
2()l060rad/s  (32,000rps),  the  sphere  was  subjected  to  a  maximum  shear  stress  about 
1 00 MPa  at  its  center.  Due  to  its  programming- friendly  concise  expression,  the  Norton 
model  was  used  to  depict  this  creep  behavior.  From  available  literature,  parameters  in 
Norton  model  at  1400  were  calculated  for  ZrB^  and  ZrB^  +  25  voi%  SiC,  which 
functioned  as  a  reference  point  for  approximating  same  parameter  at  adjacent  temperature. 
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In  order  to  gel  a  9%  equatorial  deformation  within  several  hours,  the  strain  rate  should  be 
at  the  order  of  10"'  /  sec  which  is  comparable  to  our  group’s  previous  work.  Achieving 
this  fast  creep  rate  requires  a  temperature  of  1800  for  ZrBj  and  1500  ®C  for  ZrBj  + 
Sic,  even  at  angular  velocity  as  high  as  32,000  rps.  The  simulation  results  agree  quite 
well  with  what  we  expected  and  previous  work.  The  simulations  also  show  that  this  very 
high  angular  velocity  is  necessary  to  complete  the  measurements  with  the  feasible  time 
on  the  NASA  ESL,  showing  the  importance  of  this  new  experimental  capability. 
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7.  Figure  captions: 

Fig.l .  X  direction  deformation  profile  of  XZ  plane  for  pure  ZrBj,  Generally,  material 
subjects  to  increasing  shear  force  proportional  to  its  distance  to  the  sphere  center, 

Fig.2.  Von  Mises  stress  plot  of  deformed  pure  ZrB2.  According  to  equation  {2),  Shear 
force  scales  with  square  of  its  distance  to  the  sphere  center.  This  plot  agrees  with 
theoretical  prediction. 

Fig. 3.  X  direction  deformation  profile  of  XZ  plane  for  ZrBj  +  25  vol%  SiC.  Material 
subjects  to  increasing  shear  force  proportional  to  its  distance  to  the  sphere  center, 

Fig.4.  Von  Mises  stress  plot  of  deformed  ZrB2  +  25  vol%  SiC.  According  to  equation 
(2),  Shear  force  scales  with  square  of  its  distance  to  the  sphere  center.  This  plot  agrees 
with  theoretical  prediction. 
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8.  Tables: 


Table  1,  Density  and  Young's  modulus  for  pure  ZrB^  and  ZrBj  +  25  vol%  SiC  at  room 
lemperarure.  Read  data  for  ZrBj  and  interpolate  data  for  ZrB^  +  25  voI%  SiC  from 
materials  with  close  composition  based  on  Chamberlain  ei  ai  work. 


Material 

Density 

Young's  modulus  (GPa) 

Pure  ZrBj 

6.26 

489 

ZrB2  +  25  vol,  %SiC 

5.60 

475' 

Data  with  superscript  *  means  interpolated  data. 
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Table  2,  Poisson’s  ratio  for  pure  ZrBj  and  ZrB^  +  25  voI%  SiC.  Since  Poisson’s  ratio  is 
almost  insensitive  to  additives  and  processing  conditions,  we  assume  that  Poisson’s  ratio 
is  constant  over  interested  temperature  span  of  our  simulations. 


Material 

Poisson’s  ratio 

Pure  ZrBj 

0.15 

ZrB2  +  25  voi%  SiC 

0.15 

15 


Table  3. 


Parameters  of  Norton  creep  model  for  ZrB^  and  ZrBj  +  25  vol%  SiC  as  well  as  estimated 
stress  rate. 


Material 

Temperature{®C) 

Stress  exponent 

Coef!1cent(  .s  ’ ) 

Estimated  Creep  rate  (.9*‘ ) 

HI 

1800 

3 

1x10-'' 

1.4305*10  ' 

ZrB2  +  25 

vol%  SiC 

1500 

1.3 

1x10’ 

9.1212x10  " 
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